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ABSTRACT 

We investigate the effect of galactic winds on the Lya forest in cosmological simulations 
of structure and galaxy formation. We combine high resolution N-body simulations 
of the evolution of the dark matter with a semi-analytic model for the formation and 
evolution of galaxies which includes detailed prescriptions for the long-term evolution 
of galactic winds. This model is the first to describe the evolution of outflows as a 
two-phase process (an adiabatic bubble followed by a momentum-driven shell) and to 
include metal-dependent cooling of the outflowing material. We find that the main 
statistical properties of the Lya forest, namely the flux power spectrum P{k) and the 
flux probability distribution function (PDF), are not significantly affected by winds 
and so do not significantly constrain wind models. Winds around galaxies do, however, 
produce detectable signatures in the forest, in particular, increased flux transmissivity 
inside hot bubbles, and narrow, saturated absorption lines caused by dense cooled 
shells. We find that the Lya flux transmissivity is highly enhanced near strongly wind- 
blowing galaxies, almost half of all high-redshift galaxies in our sample, in agreement 
with the results of Adelberger et al. (2005). Finally, we propose a new method to 
identify absorption lines potentially due to wind shells in the Lya forest: we calculate 
the abundance of saturated regions in spectra as a function of region width and we 
find that the number with widths smaller than about 1 A at z = 3 and 0.6 A at 2 = 2 
may be more than doubled. This should be detectable in real spectra. 
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1 INTRODUCTION 



Extensive observations of QSO spectra have provided a 
wealth of information on properties of the intergalactic 
medium (IGM). 

The Lya forest provides a measure of the power 
spectrum of density fluctuations Pm(fc) on smaller scales 
than those accessibl e by Cosmic Microwave Back ground 
(CMB) observations fepergel et ai] |2003'. 'Viel ct al.' 2004a') 
or large-scale structure surveys (e.g. .Tegmark ct al. 2004, 
IViana. Nichol fc Liddldl200l ICoUess et al.li200lh . The flux 
power spectrum P{k) of the Lya forest is directly re- 
lated to the power spectrum of density fluctuations Pm(fc), 
thanks to the existence of a direct connection between 
the Lya ab sorption and the underlying density and ve- 
locity flelds JViel et al.ll2004bl ICroft et alJbool iHui et alJ 



been estimated at red s hifts 1 < z < 4 llKim et al.ll2004L 



2001].ICroft et alJll998l) The flux power s pectrum P(k) has 

McDonald ct al.' "2004', ICroft et al.l I200I iMcDonald et alJ 

2000 , Croft et al.,,d99al for 0.001 < fc < 0.1 s km"\ corre- 
sponding to scales up to about 50 comoving Mpc. Given 
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that the data sets used by different groups were taken with 
different instruments (e.g. UVES, HIRES, LRIS) on differ- 
ent telescopes (e.g. VLT, Keck, SDSS Apache Point) and at 
different resolution and signal-to-noise ratio, the agreement 
between the results is good. 

When using the matter power spectrum Pm(fc) from 
the Lya forest to estimate cosmological parameters in con- 
junction with CMB and large-scale structure data sets, it 
is important to know precisely how P{k) is affected by 
sources of systematic errors and non-cosmological distor- 
tions. For example, the estimation of the unabsorbed con- 
tinuum level in a QSO spectrum is extremely sensitive to 
fluctuations in the cosmic UV background radiation, which 
is poorly constrained, and this affects both the amplitude 
and slope of P(k) fcidz et al . 2005, McDonald et al.ll2005l. 
IViel et alj|2004d. iHui et al.ll200L) . High column density ab- 
sorbers with damping wings and signatures of galactic winds 
may produce "local" distortions in the Lya forest which in- 
troduce deviations in the power spectrum at smaller scale s 
jFang et alJl2005l . iMcDonald eral]l200l IViel et al.ll20043) . 
The presence of metal lines blended with the Lya for- 
est affects the fl ux power spectrum at very small scales 
jKimeraOliooi)- 
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Adelberger et al. (2003) jointly analysed high resolution 
quasar spectra and spectra of Lyman break galaxies (LBG) 
near the quasar lines-of-sight. Their earliest results indi- 
cated that the Lya continuum emitted by the background 
quasars can be almost completely transmitted by the inter- 
galactic medium (IGM) near LBGs. Numerical simulations 
jDesiacaues et al.lbOo'llMaseUi et al.ll2004l.lKollmeier et alJ 
l20o3) liave been unable to reproduce this result, predicting 
instead a decrease in the Lya flux transmissivity with de- 
creasing impact parameter (the impact parameter is defined 
as the minimum distance between a galaxy and the quasar 
line-of-sight). When new results derived from a larger obser- 
vational sample were released (Adelberger et al. 2005, here- 
after A05), the improved statistics showed that the mean 
transmitted flux around LBG does indeed decrease with de- 
creasing impact parameter, as predicted by the simulations, 
but about a third of the galaxies in the full sample still show 
little or no Lya absorptio n in their proximit y . This suggests 
that the earlier result of lAdelberger et al"! (l2003ll was af- 
fected by small number statistics, although even a relatively 
small fraction of LBG's "bubbles" in absorption appears in- 
consistent with published simulations. 



2 SEMI-ANALYTIC SIMULATIONS OF 
GALACTIC WINDS 

Since we want to investigate the effects of outflows in their 
proper cosmological context, we opt for a combination of 
N-body simulations for the evolution of the dark matter 
density field and a semi-analytic model for the formation 
and evolution of galaxies. This choice gives us the possibil- 
ity to efficiently combine a high resolution in mass with a 
large simulated volume, as is necessary to study the effects 
of feedback from galaxies with a wide range in stellar mass. 
The large volume allows us to achieve a good statistics when 
investigating the effect of different galaxy populations on 
their surroundings. Our simulations assume a ACDM cos- 
mology with matter density Qrn ~ 0.3, dark energy density 
f^A = 0.7, Hubble constant h = 0.7, primordial spectral 
index n = 1 and normalisation as ~ 0.9. 

In the following subsections we briefiy describe our set 
of numerical simulations (Subsection 12.111 and our semi- 
analytic model for the physics of galactic winds (Subsection 
12.21 . However, for a more detailed description of the wind 
model, we refer the interested reader to BSW05. 



In this work, we use a semi-analytic model for the long- 
term evolution of galactic winds by Bertone, Stoehr & White 
(2005, BSW05 hereafter) to investigate the effect of galac- 
tic winds on the Lya forest. In particular, we calculate the 
flux probability distribution function (PDF) , the flux power 
spectrum P{k) and the Lya fiux transmissivity near high 
redshift galaxies, and we compare the numerical results with 
the relevant observations. We consider three different mod- 
els, a "no wind" model and two wind models corresponding 
to qualitatively different scenarios, and we try to determine 
which model better reproduces the observations. In addi- 
tion, the comparison of the "no wind" model with the two 
wind models highlights the effects of galactic winds on the 
Lya forest. To do this, we extract artificial spectr a from 
the " M3" set of high resolution N-body simulations (IStoehJ 
1200,^^ . which combines a large dynamical range with a high 
meiss resolution . The semi-analytic wind model of BSW05 
includes a new description of the dynamics of the winds as a 
two-phase process: a first phase of adiabatic evolution dur- 
ing which the gas outfiowing from galaxies is hot, and a 
second snowplough phase in which the wind material cools 
down and accumulates in a thin shell pushed by the momen- 
tum of the wind. The efficiency of winds in the simulations 
is set by a parameter (e, see Subsection 12. 211 that determines 
the amount of mass swept up by the expanding gas. 

The paper is organised as follows. In Section|5|we briefiy 
describe our set of simulations and the semi-analytic model 
for winds of BSW05. Section |3 presents our method to ex- 
tract artifical spectra from the simulations and to include 
the contributions of winds. Detailed descriptions of our cal- 
culation of the hydrogen ionisation state and of the fiux 
spectrum are given in Appendix IXI and Appendix |E| respec- 
tively. Results are presented in Section 01 for the flux PDF 
and in Section |^ for the power spectrum. In Section |S| we 
present our results for the flux transmissivity around LBG 
and we discuss the observations of A05. In Section we 
present a new method to test statistically for the signature 
of winds in spectra. Section |H| contains our conclusions. 



2.1 The M3 simulations 

We use the "MS" high- resolution N-body simulation of 
IStoehd i2003h . M3 is a resimulation at higher resolution 
of an approximately spherical region of the universe with 
average density close to the cosmic mean and diameter 52 

Mpc. The particle mass in the high resolution region is 
1.7 ■ W^h~^ Mq and the number of particles about 7 • 10^. 
The simula tions were performed using the parallel treecode 
GADGET I JSprineel et al.l2001all and 52 simulation outputs 
were stored between z — 20 and z = 0. 

The formation and evolution of galaxies is modelle d 
with the semi-analytic technique of lSpringel et al.l (l2001lJ) . 
Merging trees extracted from the simulations are used to fol- 
low the galaxy population in time, while simple prescriptions 
for gas cooling, star formation and galaxy merging model the 
processes involving the baryonic component of the galaxies. 
At 2 = 3 a total of about four hundred thousand galaxies 
are identifled and about three hundred and fifty thousand 
are present ai z — 0. About half of the galaxies at 2; = are 
field galaxies, while the rest are in groups and poor clusters. 
The two largest clusters have a total mass of about 10^* 
Mq . The recipes for the evolution of winds are implemented 
on top of this pre-existing scheme. 

2.2 The wind model 

BSW05 introduced new prescriptions for the evolution of 
galactic winds in the semi-analytic model for galaxy forma- 
tion of SpriiiEcl et al. (2001b). The most innovative feature 
of their prescriptions is a description of the dynamics of out- 
flows as a two-phase process: a pressure-driven adiabatic 
expansion followed by a momentum-driven snowplough. 

Previous simulations of galactic winds in a cosmo- 
logical context usually assume that winds are either 
pre ssure-driven or mom e ntum- driven. The wind model 
of ISpringel fc Hernauisd (l2003l) mimics the behaviour 
of an adiabatic bubble expanding into the surround- 
ing medium. However, since metal cooling is neglected. 
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ISpringel fc Hernauisd i2003l) predict too high temperatures 
for the ejected metals and they are unable to reproduce the 
observe d C IV ab sorption in QSO spectra ijAgxiirre et alJ 
boos) . Theuns et al . (2002) find that winds have little ef- 
fect on the H I absorption statistics, because the hot bubbles 
mostly escape into the voids. Semi-analytic models of galaxy 
formation that include winds us ually describe outfl ows as 
momentum-driven snowploughs f A guirre et alj|200iri . 

We do not resolve the first phases of the wind evolution, 
when a superbubble sweeps through the ISM of a galaxy. 
Instead, we follow the long-term evolution of winds once 
they have escaped the visible regions of galaxies. We make 
the simplifying assumption of spherical symmetry for the 
wind evolution. This may seem a rough approximation at 
2 ~ 0, where outflows are mostly bipolar, but it is a good 
assumption at higher redshifts, where outflows seem to have 
a mo re spherical geometry ( Shaplcv ct al. 2003, Rupkc ct aL 
|2003). 

Galactic winds are modelled as uniform pressure-driven 
bubbles of hot gas emerging from star-forming galaxies. The 
assumption that winds are adiabatic at blow-out is mo- 
tivated by o bservations of galaxies i n the local universe. 
For ex ample. iHoopes et all 1*2003) and lStrickland fc Steven^ 
l|2000fl observe no energy losses in M82 through radiative 
cooling of the coronal (T ~ 10^'^ K) and the hot (T ~ 10^ K) 
phases of the wind, supporting the idea that the early evo- 
lution of this wind is nearly adiabatic. The adiabatic phase 
of the wind evolution is described by the equation for the 
conservation of energy E JOstriker fc McKeelll988h : 



at 2 

1 2 GA4 



VoPo } , (1 



where R and Vs are the radius and the velocity of the shock, 
Mw and the mass outflow rate and the outflow velocity 
of the wind, po, Po and «□ the density, the pressure and the 
outward velocity of the surrounding medium and Mh the 
total mass internal to the shock radius. The entrainment 
fraction parameter e defines the fraction of mass that the 
wind sweeps up while crossing the ambient medium. Here, 
we call "entrained" gas the ambient gas which has mixed 
into the hot bubble phase either through turbulent mixing of 
shocked diffuse ambient gas or by evaporation and ablation 
of the filament gas, most of which (1 — e) continues falling 
onto the galaxy. The latter process is similar to the loading 
of ISM mass onto stellar winds and supernova blastwaves. 

The adiabatic phase is terminated when the loss of en- 
ergy by radiation becomes substantial and most of the en- 
ergy transferred to the swept-up gas is radiated away. When 
the cooling time of the hot bubble becomes shorter than the 
age of the wind, a thin shell of cooled gas forms near the bub- 
ble's outer boundary and continues to expand pushed by the 
momentum of the wind. The snowplough phase is described 
by the equation for the conservation of momentum: 



GMh 



eAttR^ [Po + PoVo (Ws — Vo)] , 



(2) 



where m is the mass of the shell. 

In our simulations both bursts of star formation and 
quiescent star formation can power winds, since we do not 



explicitly restrict the star formation rate required for a 
galaxy to blow a wind. It is not possible to predict a priori 
when a wind will escape the gravitational pull of a galaxy, 
since its evolution and its final fate are linked to several fac- 
tors, like the star formation and the mass accretion history 
of the galaxy, the potential well of the dark matter halo in 
which it expands, the amount of mass accreted from the 
wind and from the IGM and so on. 

In general, a wind receives energy from the starburst 
and is slowed down by the gravitational attraction of the 
central galaxy and by the ram pressure of the ambient 
medium. Thermal pressure effects are included consistently 
inside cluster haloes, but are neglected in the IGM. If the 
entrained mass is small (e.g. e ^ 0.1), a large fraction of 
the bubble or shell mass consists of supernova ejecta and 
shocked ISM, which flow out from the galaxy with a veloc- 
ity often much larger than the escape velocity of the halo. 
Since little energy or momentum has to be spent by the 
wind to accelerate the entrained mass, the shock velocity is 
less sensitive to energy losses by pressure and gravity. Such 
a wind is thus more likely to break free from the halo than 
more mass-loaded winds. When mass-loading is substantial 
(e.g. e = 0.3), a significant part of the wind energy is con- 
sumed to accelerate the entrained gas and the expansion 
slows down considerably. If the amount of swept-up mass is 
large compared to the mass initially ejected, the shock ve- 
locity may drop below the escape velocity from the galaxy 
and the wind collapses back onto the galaxy. 

At high redshift winds tend to be mostly momentum- 
driven, while at lower redshifts bubbles are much more likely 
to remain adiabatic. This is partly due to the higher mean 
density of the Universe at high z and partly to a lower en- 
ergy input from star formation, which results in lower bub- 
ble temperatures and shorter cooling times immediately af- 
ter blowout. The transition from pressure-driven bubbles to 
momentum-driven shells may be a first hint that a wind is 
not powerful enough to escape the galaxy's attraction. In 
fact, pressure-driven winds are overall more likely to escape 
galaxies than momentum-driven ones. 



3 SIMULATING THE Lya FOREST 

In this section we show how we extract sets of artificial spec- 
tra from our simulated region. Since we use a semi-analytic 
model associated to N-body simulations, we do not have a 
priori all the information about the gas properties provided 
by a gas-dynamical simulation. We therefore need to esti- 
mate quantities like the gas temperature and density, which 
are necessary to build the spectra, by applying approximate 
prescriptions. We show how in Subsection l3.1l In Subsection 
13.21 we describe how we include the effects of winds in the 
calculation of the spectra. 

We extract lines-of-sight (LOS, hereafter) from our 
simulated re g ion following the numerical scheme of 
iTheuns et alJ 1^9^. We describe our assumptions about 
the ionisation state of the gas along the LOS in Appendix 
IXI In Appendix^ we summarize our prescriptions for the 
integration of the relevant quantities along the LOS and for 
the calculation of the optical depth and the flux spectrum. 
Each spectrum i s normalized to the mean observed flux of 
iKim et all i2002l) . Gaussian noise is added to the simulated 
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spectra assuming a fixed signal-to-noise ratio of 50, typical 
of HIRES spectra. Gaussian read-out noise with variance 
0.004^ is also added to each pixel. 

We extract synthetic spectra from three different sets 
of simulations: 1) a "no wind" model, in which we do not 
include the effects of galactic winds; 2) a wind model with 
low mass loading efficiency and entrainment fraction param- 
eter e = 0.1 (hereafter "eOl"); 3) a wind model with high 
mass loading efficiency and entrainment fraction parameter 
e = 0.3 (hereafter "e03"). For each set of simulations we 
extract spectra at 2 ~ 3 and z ~ 2. 

For each wind model, we then extract two sets of spectra 
from our simulated region: i) one set along random directions 
and ii) one set imposing that the LOS intercept specially se- 
lected galaxies at random distances in the interval (0, 10) 

Mpc comoving. We use the first sample of spectra to 
calculate the fiux probability distribution function (PDF) 
and the power spectrum P{k) and the second to investigate 
the effects of winds on the surroundings of galaxies. The sec- 
ond set of spectra was selected to match the observations of 
A05 as accurately as possible. In table 3 of AOS are listed the 
properties of the galaxies observed by NIRSPEC with impact 
parameters (e.g. distance galaxy-LOS) smaller than 1 
Mpc comoving. According to these data, we select galaxies 
in our simulation which have: 1) stellar masses larger than 
10^ Mq; 2) SFR larger than 1 Mq/yv: these two values 
are the minimum stellar mass and the minimum SFR esti- 
mated by A05. In total we extract a sample of about 4400 
LOS at 2 ~ 3 and about 3300 LOS at 2 ~ 2: one LOS for 
every galaxy with the desired properties. We make no re- 
quirements about the presence of winds, but we do require 
that both the selected galaxy and the LOS fall inside our 
high resolution region. 

The set of random spectra has been selected assuming 
that each LOS passes at a distance of less than half the 
radius of the high resolution region from the centre of mass 
of the simulation Xc, that is d (LOS, Xc) < 13 Mpc. This 
is because we want to extract simulated spectra of minimum 
length 45 Mpc comoving. 

3.1 Density and temperature of the gas 

The gas density in our simulated region is assumed to follow 
the dark matter density distribution. This is not a very satis- 
factory description for large matter overdensities, like galax- 
ies or galaxy groups, but it is sufficiently accurate for the gas 
in low density regions, which represents about 70% of the 
mass at 2 ~ 3 and which is responsible for most of the Lya 
forest. The density of the dark matter pdm is calculated by 
SPH interpolation over 26 neighbours and the density of the 
gas pgas is recovered by assuming that pgas = .fbaryon • Pdm, 
with /baryon the global baryon fraction. 

The temperature of the "gas" particles is calculated 
separately for high density particles bound to haloes and 
low density unbound particles. This double treatment is 
necessary since the physical conditions of the gas in the 
two cases are different. The intergalactic gas is mostly un- 
shocked, cools adiabatically and is photoionised by the UV 
background radiation. On the other hand, the gas in haloes is 
shock-heated because of the structure formation process and 
its temperature is determined by the gravitational potential 
of the dark matter halo in which it resides. We therefore 
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Figure 1. The temperature of the gas in M3 at 2 = 3. Gas 
that obeys the equation of state in eq.|3] piles up on the straight 
Hne, while the shocked gas in haloes is spread out in a cloud-like 
region. 



divide our particles in two subgroups, the "halo" particles 
and the "IGM" particles, and calculate their temperatures 
according to different prescriptions. 

For "halo" particles we mean all those particles which 
belong to bound structures and we assume that the temper- 
ature of_th£S£_2ajticles_emi^^ the virial temperature of the 
halo JWhite fc Frenklll99 j) : 



Tvir — 



2kB 



35.9 



T4 



km s 



K, 



(3) 



where Vc is the circular velocity of the halo, p is the mean 
molecular weight of the gas, mp is the proton mass, pmp 
the mean particle mass and fee the Boltzmann constant. 

The low density "IGM" particles represent the dif- 
fuse intergalactic medium and usually have overdensities 
S < 10 — 100. For these particles, we calculate the temper- 
ature from the equation of state of jHui fc Gnedin. L1997t) . 
which approximates a power law: 



To (1 + S) 



7-1 



(4) 



where 7 — 1 ~ 1/1.7. To is the temperature of the IGM 
at the mean density and can be recovered under the as- 
sumption that the gas is in photoioni sation equilibrium 
jHui fc Gnedinl '1997'. 'Schave et al."2000'). In this work, we 
use the results of Bolton ot al. ( 2004) to define To as a func- 
tion of redshift. 

In figure^we show the temperature of the gas in M3 at 
2 = 3. The low density gas obeying the equation of state ||1J 
is distributed along a straight line, while the high temper- 
ature gas in haloes is spread into a cloud-like region. The 
horizontal lines crossing the cloud are due to gas with dif- 
ferent densities residing in the same halo, which is assumed 
to behave as an isothermal gas sphere. Our prescriptions for 
the gas temperature are unable to reproduce the cold neu- 
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tral gas in galaxies and galaxy haloes. This will partly affect 
our results and we will show how in the following Sections. 

3.2 The matter displaced by winds 

To include the effects of winds in the calculation of a spec- 
trum, we proceed in the following way. As seen in Section 
m\ winds can be described as bubbles of hot gas or as thin 
shells of cold material accumulating at the edge of a cooled 
cavity. These two phases of the wind evolution produce dif- 
ferent conditions for the material displaced by winds and a 
separate treatment of bubbles and shells is therefore needed. 

Since winds modify the spatial distribution of the gas 
near galaxies, we need to identify how much mass is dis- 
placed, which particles can be associated with the displaced 
matter and where this matter moves. All the particles resid- 
ing outside winds remain unaffected. To do this, we identify 
all the gas particles positioned inside winds. We then calcu- 
late the total amount of mass Min inside winds as predicted 
by the semi-analytic model: this is equal to the sum of the 
stellar and cold gas mass of galaxies plus the halo or IGM 
mass inside the wind radius. The total mass contained in 
bubbles, shells and cavities is Mout = Yl (-^w + nic), with 
Mw the mass of the wind ejecta and me the swept-up mass 
in each wind. We finally define F = Mout /Min as the frac- 
tion of the baryonic mass in the wind-affected region that 
belongs to bubbles and shells. 

Once F is calculated, we determine which individ- 
ual particles are "wind" particles and which ones are not; 
"wind" particles are ignored when integrating quantities 
along the LOS and alternative prescriptions are used to cal- 
culate the optical depth contributed by winds, as we will 
show in the following. Since there is no exact or unique way 
to identify the "wind" particles, we decided to flag as "wind" 
a fraction F of all particles in the wind-affected regions, 
speciflcally those with the lowest estimated densities. This 
approximation implies that the particles we remove from the 
list are generally residing far from the centre of haloes, while 
the particles representing the galaxy and the unperturbed 
ambient medium tend to cluster towards the centres. How- 
ever, since particles have a finite volume, defined by their 
SPH smoothing radius hi, the lowest density particles are 
also the ones spread over the largest volumes, and so most 
likely to contribute to the LOS. Contributions from matter 
"in" galaxies is mostly unaffected. 

In our calculations, we first determine the density and 
the density-weighted temperature and velocity contributed 
by non-wind particles to each bin j along a LOS (Appendix 
[Hi . We then add in the contribution of the winds themselves. 
Finally, we calculate the optical depth along the LOS and 
the fiux spectrum. In Subsections 13.2. H and l?!2.2l we show 
how we estimate the contributions of bubbles and shells re- 
spectively and we discuss a few examples of synthetic spec- 
tra. 

3.2.1 Contribution by bubbles 

To calculate the contribution of bubbles to the optical depth 
along a LOS, we assume that the material outflowing from 
the galaxies and the shocked ambient medium are efficiently 
mixed and uniformly distributed inside the bubble. This as- 
sumption of a uniform bubble is reasonable (although not 
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Figure 2. Example of spectrum crossing a pressure-driven bub- 
ble. The blue line shows the unperturbed "no wind" model, the 
red line the e03 wind model. From top to bottom, the different 
panels reproduce the flux spectrum, the optical depth t of the 
gas along the line of sight, the fraction of neutral hydrogen Xn i, 
the total density p, the temperature T and the peculiar velocity 

^pcc • 

exact), because the sound crossing time of the gas is shorter 
than the dynamical time of the wind. 

The semi-analytic description of the wind evolution 
gives information about the mass, radius, velocity and tem- 
perature of the wind. The temperature of the bubble Tb is 
determined by the energy balance in the wind (cfr. eq. 4 of 
BSW05). Most temperatures are of order 10^ to few times 
10^ K, which means that the bubble material is coUisionally 
ionised. We calculate the fraction of neutral hydrogen Xh i 
in a bubble by assuming coUisional ionisation equilibrium 
( Appendix IXt. 

The neutral hydrogen density and the density-weighted 
temperature and velocity that a bubble contributes to each 
bin j which intercepts the bubble, are: 

pa = "^-'^H iPh, (5) 
(pr)H J (j) = a'Xn iPbTb, (6) 

(P")h I (j) = «^^H I/SbWpocb, (7) 

where a is the scale factor, pb the bubble density and Vpac 
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the projection of the pecuhar velocity of the wind on the 
LOS. 

In Fig. |21we show an example of the build-up of a spec- 
trum along a LOS that crosses a bubble. From top to bot- 
tom, the panels show the flux spectrum, the optical depth 
T, the fraction of neutral hydrogen Xh i, the total density 
p, the temperature T and the peculiar velocity Vpcc of the 
gas along the LOS. The blue line represents the unperturbed 
quantities calculated from the "no wind" model, while the 
red line is for the perturbed quantities from the e03 model. 
In this spectrum, the LOS intersects a large bubble at a 
small distance from the source galaxy, but does not intercept 
the galaxy itself. The ionised gas inside the bubble produces 
a decrease of the optical depth, due to its high temperature 
and ionisation state. The "Z" shape in the velocity distribu- 
tion of the gas along the LOS (bottom panel) is due to the 
outflowing wind material, Vs ~ 200 km s~^, which modifies 
the velocity field around the galaxy. 

The flux spectra and the optical depth differ slightly in 
the two cases because of normalization effects: this happens 
when winds modify the mean optical depth along the LOS 
by more than a few percent and an iteration or two is needed 
during the normalization process. However, the final differ- 
ence between the "wind" and the "no wind" models after 
normalization is always less then a few percent, even in the 
most extreme cases. 



3.2.2 Contribution by shells and cavities 

The effect of shells on the opacity of the gas along a LOS 
is substantially different from that of bubbles. Both the 
shell and the cavity evacuated in its interior are cold, and 
have a more complex density structure. When winds cool 
down, part of the outflowing wind mass and all the swept- 
up mass accumulate in a thin shell. A cavity is left behind, 
in which the low-density, cooled wind material that has not 
yet reached the shell flows outwards at the wind velocity. 
The mass that accumulates in the shell occupies a small 
volume and its density can become significantly higher than 
the density of the surrounding medium. 

We make the assumption that the temperature of shells 
and cavities is constant and equal to Ts — 10* K. This as- 
sumption is motivated by the fact that, after the wind cools 
down and all the energy is radiated away, the outflowing 
material reaches photoionisation equilibrium with the UV 
background and will maintain this temperature as long as 
the density of the shell is not too high. 

Our semi-analytic model does not give an explicit de- 
scription of the density structure of shells and cavities, but 
we can reconstruct it by solving the appropriate set of shock 
conditions. We use the Rankine-Hugoniot jump conditions 
for a non-relativistic shock to calculate the pressure Ps of 
the gas in the shell. If we assume that the velocity of the 
gas behind the shock is equal to the shock velocity Vs, we 
obtain: 



(8) 



where Pigm, Vigm and -Rgm are the density, the velocity and 
the pressure (thermal plus ram) of the ambient medium. 
The density of the shell ps is then given by the equation of 
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Figure 3. Example of spectrum crossing a momentum-driven 
shell. The quantities shown are the same as for Fig. |21 for the 
unperturbed "no wind" model (blue line) and for the e03 wind 
model (red line). 



state: 



p.mpPs 



(9) 



The radius Rc of the cavity internal to the shell is given by: 



itr — it 



3m 
47rps 



(10) 



and the thickness of the shell is Ra = R — Rc, with m the 
shell mass. The density of the cavity is determined by the 
wind mass which has not yet reached the shell Mc and is 
given by pc = 3Mc/4t7R^. 

We estimate the contribution of shells and cavities to 
the neutral hydrogen density and the density-weighted tem- 
perature and velocity along the LOS by using eqs. ^ 
and O , where the index "6" has been changed with "c" and 
"s" . Shells affect only a few bins on the outer edge of a wind 
penetrated by a LOS. The number of the bins is determined 
by the thickness of the shell, which is normally as large as 
the linear dimension of few bins (typically, the bin width is 
about 20-26 kpc). The remaining bins inside the shell 
get contributions from the cavity gas. 

In Fig.|3we give an example of the build-up of a spec- 
trum that crosses a shell. The quantities shown are the same 
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Figure 4. Flux PDF for the Lya forest at z = 3. The numeri- 
cal re sults are compared to the observations of lMcDonald et alj 
J200(l) . 

as in Fig. |5| The shell crossed by the LOS has a density 
about ten times the density of the ambient medium, but the 
fraction of H i in the shell is about 100 times higher. The 
cold neutral gas in the shell produces two almost symmetric 
spikes in the optical depth, which appear as two narrow and 
saturated absorption lines in the flux spectrum. The tem- 
perature of 10* K that we assume for the shell determines 
the small width of the lines, while the high neutral fraction 
produces the saturation. 

The LOS in Fig. |21 and |3 show examples of the effects 
of winds in spectra. In general, a line of sight may not cross 
any wind at all or it may cross several, in which case the 
spectrum may become more confused. 



4 FLUX PDF 

We calculate the flux probability distribution function 
(PDF) as a first comparison between our simulated spec- 
tra and the observations. We use about 600 spectra from our 
random sample at z = 3 and we compare the results with the 
observations of McDonald ot al. (2000). Fig.0]shows our re- 
sults for the "no wind" model and for the eOl and e03 wind 
models. The agreement between the observations and our 
three sets of simulations is in general very good and devia- 
tions from the observed values are significant only at fiuxes 
F ~ and F > 0.9. However, the observed values corre- 
sponding to F ^ 0.9 are not fuUy reliable, because they can 
be strongly affected by continuum fitting errors. 

In general, the e03 model gives almost exactly the same 
result as the "no wind" model, while a small but more sig- 
nificant difference is evident with the eOl model. This may 
be an indication that the weak winds predicted by the e03 
model are unable to affect the fiux PDF of the Lya forest, 
while more efficient winds as in the eOl model may slightly 
modify it. However, in all cases the difference between the 



Figure 5. The flux power spectrum computed from synthetic 
spectra at z ~ 3. The numerica l results are compared to the 
observ ations of iKim et al 1 i2004) at z ~ 2 54, iMcDonald et alj 
l2000l) at 2 = 3 and iMcDonald et all l2004) at 2 = 3. The syn- 
thetic spectra have been calculated using the samples of randomly 
selected LOS for the "no wind", eOl and e03 models. 

numerical results and the observations are tiny, and not suf- 
ficient to rule out any of the feedback models. 



5 POWER SPECTRUM 

To calculate the po wer spectrum of t he flux, we use the esti- 
mator proposed bv lHui et al.l ll200 j) . Since the power spec- 
trum is sensitive to changes in the mean s pectr al flux {e~^), 
which is a function of redshift. lHui et al] ||200J) suggest cal- 
culating the power spectrum of the quantity: 



where r is the optical depth and the amplitude of the 
fluctuations. Fig.l^shows our results for the power spectrum 
of the transmitted flux P{k), calculated in this way from the 
sets of random spectra at z = 3 for the "no wind", eOl and 
e03 models. All the spectra for the three wind models have 
been calculated assuming the same relationship between the 
gas density and temperature, as given in Subsection 13.11 
This is important, because the flux spectrum is sensitive to 
changes in the density-temperature relationship: the same 
gas conditions ensure that our results show variations due 
to winds and not to simulation artefacts. 

The nu merical res ults are compared to the obser- 

vatioi is of iKim et all ||2004) . IMcDonald et al.l (l2000h 
and 'McDon aM et alJ 11200411 . iKim et alj (12004) and 
IMcDo nald ct al] ll200Cl) compute the power spectrum from 
few high resolution quasar spectra taken with the VLT and 
Keck telescopes respectively, while IMcDonald et al] (l2004h 
use a sample of about 3000 medium resolution spectra from 
the Sloan Digital Sky Survey (SDSS). The lower resolution 
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of the SDSS spectra limits the estimate of P{k) in the range 
0.0013 — 0.04 s km"'^, while higher resolution spectra allow 
to estimate P{k) in the broader range 0.001 — 0.2 s km~^. 
At fe > 0.3 — 0.4 s km~^ P{k) is dominated by noise and 
flattens out. 

In Fig. |S] there is agreement between the "no wind" 
model and the eOl and e03 wind models at fc < 0.1 s km~^, 
with deviations appearing only at smaller scales. At A; > 0.1 
s km~^ the three numerical models diverge and, broadly 
speaking, the higher the wind efficiency, the more power 
is contributed at smaller scales. This is consistent with our 
previous results for the wind filling factor and the fraction of 
IGM mass in winds (BSW05), which predict that, although 
winds may not heavily modify the statistical properties of 
the Lya forest, the largest deviations should be found in the 
eOl model. In addition, the fact that more efficient winds 
contribute more power on small scales is consistent with 
the presence of the narrow and saturated absorption lines 
created by shells, which are more numerous in the eOl model. 
These lines contribute power to scales smaller than about a 
hundred kiloparsecs. 

Unfortunately, the spectral region corresponding to k > 
0.1 s km^^ is also the region where unidentified metal 
lines (principally O VI and C IV) blended with the Lya 
forest might give a substantial contribution to the power. 
iKim et all (l2004l) estimate that between 20 and 70% of the 
power is contributed by metal lines for k > 0.1 s km~^. Al- 
though we are tempted to suggest that both our wind models 
seem to fit the data better than the "no wind" model, we 
cannot definitely rule out the "no wind" model itself, since 
our current Lya spectra do not include metal lines and we 
are unable to estimate well how much power is contributed 
by them. 

The good agreement between the two wind models and 
the "no wind" model at A: < 0.1 s km~^ is an indication 
that the effects of galaxy feedback on the Lya forest, and in 
particular of galactic winds, should not affect the estimation 
of cosmogical parameters. 



6 FLUX TRANSMISSIVITY AROUND LEG 

In this Section we present our results for the mean Lya flux 
transmissivity near Lyman break galaxies and we investigate 
the effect of winds on the absorption spectra at small impact 
parameter. We use our sample of galaxy-selected spectra at 
z = 3 and z — 2 for the "no wind", eOl and e03 models 
(Section 

Fig.|S|shows our results for the mean transmitted flux as 
a function of the impact parameter R. The numerical results 
are compared to the observations of AOS. The analysis of our 
full set of spectra conflrms the observation that the mean 
Lya transmissivity near LBGs decreases for small impact 
parameters. Our spectra yield a smaller decrease than seen 
in the observations: this is partly due to the fact that our 
semi-analytic prescriptions are unable to describe the cold 
neutral gas in haloes and simply assume that this gas is 
at the virial temperature and so highly ionised. In addition, 
the mean flux (F) = 0.792 in the simulated spectra is higher 
than in the data, where (F) — 0.765. This is not a severe 
problem, however, because the difference is small and all our 
models do predict a non-negligible decrease in the Lya flux. 
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Figure 6. The mean flux transmissivity around LBG as a func- 
tion of the impact parameter. Results are presented for a sample 
of about 7700 simulated galaxies with random impact parameters 
(lines) and for subsamples of galaxies blowing winds with radii 
larger than the impact paramater (diamonds). The observational 
data are from A05. 
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Figure 7. The total number of galaxies and the number of wind- 
blowing galaxies with wind radii larger than the impact parameter 
as a function of impact parameter in the simulated sample. 



A more detailed analysis shows that the Lya flux is al- 
most completely transmitted around a significant fraction of 
the galaxies in our model and that the gas in these regions 
has an extremely low optical depth. This is shown in Fig. |^ 
by the diamonds in the upper left corner, which represent the 
mean transmitted flux along those LOS that actually inter- 
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Figure 8. Two sets of parallel lines of sight intercepting isolated field galaxies blowing winds into the IGM at impact paramaters 
100 < b < 900 comoving kpc. The left panel shows how the absorption from neutral gas surrounding a galaxy is quenched inside a 
hot pressure— driven bubble. The right panel shows the effect of a momentum-driven wind and the absorption by the cold gas in the shell. 
The red line represents the spectrum perturbed by winds from the e03 model, while the blue line the unperturbed spectrum from the 
"no wind" model. The galaxy in the left panel has = 1.6 ■ 10^ Mq and is blowing a bubble with radius R = 257 h^^ (physical) 
kpc, velocity 274 km s~^ and temperature T = 2.2 ■ 10^ K. The galaxy in the right panel has Af* = 9.2 • 10® Mq and is pushing a 
shell with radius _R = 272 (physical) kpc and velocity 103 km s~^. 



cept wind bubbles. These galaxies show transmitted fluxes 
up to 50-60% higher than the mean flux found for the fuU 
sample. 

In Fig.|7|we show the total number of LOS and the num- 
ber of intercepted wind bubbles in the eOl and e03 models 
as a function of impact parameter. Since we assign random 
impact parameters when extracting the LOS, the distribu- 
tion of the total number of galaxies is close to uniform. The 
eOl model produces wind-blowing galaxies with bubble radii 
as large as 2.5 Mpc, while fewer winds and maximum 
radii of 1.5 Mpc are produced by the e03 model. The 
fraction of galaxies for which the LOS intercepts a bubble 
for impact parameters below 1 Mpc is about 7% for the 
e03 model and about 40% for the eOl model. These numbers 
are in good agreement with the finding of A05 that about a 
third of the galaxies in their sample show almost complete 
transmission of the quasar Lya flux at the LBG redshift. 

Winds either evacuate cavities or drive bubbles of 
hot gas into their surroundings. Although the physics of 
pressure-driven and momentum-driven winds is different, 
in either case the optical depth of the gas may be reduced 
inside the wind, and the flux increased. In the case of adi- 
abatic winds, the hot gas in bubbles is highly ionised, bub- 
ble temperatures being normally higher than 10^ K. In case 
of momentum-driven winds, the swept-up gas accumulates 



into a thin shell, while a cavity is evacuated in the interior, 
where only the cooled wind is flowing outwards. The density 
of the residual gas inside the cavity is tiny and its contribu- 
tion to the optical depth along the LOS is generally negli- 
gible. Occasionally another absorber external to the galaxy 
may intervene to produce an uncorrelated absorption line. 

In Fig.lHlwe show the effects of winds on artificial spec- 
tra at various impact parameters from source galaxies in the 
e03 wind model. For clarity, we have chosen two isolated field 
galaxies with wind radii larger than 1 comoving Mpc 
and with no other wind-blowing galaxies in their proximity. 
In regions where galaxies are more clustered, winds overlap 
and their signatures in spectra mix and become less easily 
identifiable. The left panel shows the effect of a pressure- 
driven bubble of hot gas that ionises the gas surrounding 
the source galaxy (see also Fig. I^J. Since there are no inter- 
vening absorbers along the LOS, the bubble increases the 
fiux transmissivity in the spectral region corresponding to 
the galaxy position. The right panel shows a momentum- 
driven shell that evacuates a cavity and accumulates mate- 
rial at its edges (see also Fig. |^. The symmetric pattern of 
two narrow and saturated absorption lines is clearly visible. 
The separation between the two absorption lines and the 
pixel corresponding to the minimum distance to the galaxy 
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Figure 9. The distribution of galaxies with impact parameter 
6 < 1 Mpc as a function of the flux decrement observed in 
the spectra. The observational data are from AOS. 



is roughly equal to the velocity of the shell projected onto 
the LOS. 

In Fig. |5]we show the distribution of galaxies with im- 
pact parameter smaller than 1 Mpc as a function of the 
flux decrement. The numerical results are compared with 
the data of A05 and normalized to the total number of ob- 
jects in the observational sample (e.g. 24 galaxies). The flux 
decrement is defined as Fd = 1 — F. 

AOS claim that numerical simulations are unable to re- 
produce the high number of galaxies with _Pd ~ found in 
their data. In Fig.|5]we show that this is not the case for our 
model, independent of the wind parameters. In fact, the dis- 
tribution of galaxies as a function of the flux decrement is an 
alternative way to represent the PDF of the flux in the spec- 
tra and, if the impact parameters of the LOS are randomly 
chosen in the interval (0, l)h~^ Mpc, then the distribution 
of galaxies in Fig. |U] should qualitatively look like the PDF 
of the flux in Fig. 2] where the 3;-axis has been inverted 
from F to Fu = 1 — F. The simulation result presented by 
A05 that almost no galaxies have Fo ~ could be under- 
stood if the impact parameter of the galaxies in the sample 
is close to or smaller than the virial radius. In this case, a 
LOS would intercept the reservoir of neutral hydrogen in the 
halo/disk of the galaxy and a large flux decrement should 
be expected. Although our galaxy distribution explains the 
relatively high number of galaxies with Fu ~ 0, a more de- 
tailed comparison between our data and the observations is 
inconclusive. In fact, the 24 observed galaxies are too few to 
deflne a distribution calculated over 10 bins. For compari- 
son, we have almost 800 galaxies in our model sample. The 
eOl and e03 wind models predict a slightly larger number 
of galaxies with small flux decrements than the "no wind" 
model, in agreement with what we showed in Fig. |S| 



7 GAP STATISTICS 

We showed in the previous sections that winds do leave spe- 
ciflc signatures in the Lya forest, although they may not 
strongly affect the flux PDF and the flux power spectrum. In 
this section, we present a new method to test statistically for 
the presence of wind absorption lines in the Lya forest. The 
method is based on two flndings: first, when a LOS crosses 
a shell, a narrow and saturated absorption line is created; 
secondly, these absorption lines due to shells appear to have 
roughly all the same line width and steep profile. We look 
for shell absorption lines in our artifical spectra by identify- 
ing all the spectral regions where absorption is complete. To 
select such regions, we require that the flux is smaller than 
F < 0.1 throughout. We estimate the width in wavelength 
of each region, which we call Ds, and finally we calculate 
the distribution of the saturated regions as a function of the 
region width Ds- 

Our results are shown in Fig. 1101 for z = 3 and z = 2. 
The distribution of saturated regions is not affected by winds 
except for very small values of Ds, where absorption lines 
due to shells in the wind models almost double the number 
of saturated regions with 0.3 < Ds < 1 A at z = 3. At 
z = 2, where the number and the width of saturated regions 
shrinks significantly because of the cosmic expansion and 
the increased intensity of the UV background radiation, the 
wind models produce almost three times as many lines with 
0.1 < -Ds < 0.6 A as the "no wind" model. This analysis 
cannot be directly translated into a distribution of lines as a 
function of column density, because the width of saturated 
regions Ds is by no means a direct estimate of the column 
density of the absorber. However, our results show that the 
presence of absorption lines due to shells should significantly 
increase the number of narrow saturated absorption lines. 

We then tried a similar analysis to test for the presence 
of bubbles. In analogy to what we did for shells, we start 
from the result that bubbles may enhance the transmissiv- 
ity of the Lya fiux in their interior, sometimes completely 
washing out the absorption by neutral gas. We therefore look 
for "gaps" in the Lya forest and we identify all the regions 
with fluxes F > 0.9, calling their width Dg. We then calcu- 
late the cumulative distribution of the gaps as a function of 
this width. If bubbles had a large effect on the distribution, 
we would expect to find an increase in the number of gaps 
at the largest widths. The effect should be stronger at 2 = 3 
than at 2 = 2, since the gap width in the absence of bub- 
bles tends to be larger at lower redshift because of the lower 
mean opacity of the IGM. 

As can be seen in Fig. 1111 we do find such an effect, 
which is particularly significant for the eOl model. This 
model predicts up to 50% more gaps with 10 < Dg < 20 
A and about 5-20% more gaps with 5 < Dg < 10 A at both 
redshifts. For 10 < Dg < 20 A the e03 model shows a max- 
imum deviation of about 20% from the "no wind" model 
at z = 3 and almost no deviation at 2 = 2. At smaller gap 
widths there is no difference from the "no wind" model. The 
larger deviations from the "no wind" model in the results for 
the eOl model can be understood when two main properties 
of the winds are considered: i) the total number of wind- 
blowing galaxies is about ten times larger in the eOl model 
than in the e03 model; ii) the fraction of pressure-driven 
bubbles in the eOl model is two times larger than in the e03 
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Figure 10. The number distribution of saturated regions as a function of region width. Results are shown at 2: = 3 (left panel) and 
2 = 2 (right panel). Shells contribute a large number of saturated lines with widths smaller than about 1 A. 
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Figure 11. The cumulative distribution of gaps in spectra as a function of gap width. Results arc shown at z = 3 (loft panel) and z = 2 
(right panel). The distribution of gaps in the eOl model is mostly affected for gap widths 10 < Dg < 20 A while the e03 model shows a 
much smaller deviation from the "no wind" model. 



model. The maximum radius of bubbles in the eOl model is 
also larger than in the other model. 

The effect of bubbles may be difficult to distinguish 
in real spectra. The unperturbed cumulative distribution 
of gaps is unknown and no direct comparison is available. 
Metal lines may somewhat alter the distribution itself, but 
this problem could in principle be overcome, because only 
high column density metal lines not blended in the Lya for- 
est would affect the distribution, and these should be rela- 
tively straightforward to identify and remove. Althought the 



first problem is more difficult to address, we believe that 
the effect of bubbles in spectra could still be successfully 
investigated. In fact, in Fig. Illl the cumulative distribution 
of gaps in unperturbed spectra closely follows a power law 
relation, while wind models introduce a "bump" in the dis- 
tribution for 10 < Dg < 20 A. This deviation from a power 
law, as we mentioned in the previous paragraph, can be as 
high as 50% for the eOl model and it could be used to test 
for the presence of winds in spectra: on the basis of our re- 
sult, if the cumulative distribution follows a power law, most 
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likely there are no bubble signatures in spectra; on the other 
hand, if a bump does appear in the cumulative distribution 
at 10 < Dg < 20 A, then hot winds may leave a detectable 
footprint in the forest. 



8 CONCLUSIONS 

In this paper, we have used immerical simulations to inves- 
tigate the effect of galactic winds on the Lya forest. Our 
main conclusions can be summarized as follows: 

(i) The flux PDF calculated from simulated s pectra 
agrees well with the flux PDF of iMcDonald et alJ (l2000l) 
at z = 3 for all three wind models. The effects of the winds 
are barely discernible the PDF. The small deviations ob- 
served between the numerical results and the observational 
data are insufficient to rule out any of the models. 

(ii) Our simulated power spectra of the L ya flux at 
z = 3 agree well with the obser vatio nal results of [Kim e^ 
l|2004 .lMcDonal d et alJ J2000l) and lMcDonald et all i2004l'l 
on all scales, with t he eOl wind model closely matching 
the iKim et ah) l|2y(|^ 4|l data and the e03 model best fitting 
the data of lMcDonald et aP i2000l) on small scales. The 
"no wind" model tends to underestimate the power on the 
smallest scales. Winds sensibly increase the power on scales 
k > 0.1 s km~^ at 2: = 3, when both the volume filling fac- 
tor /v and the fraction of IGM mass in winds /m are small. 
The overall effect is not large, however, and deviations af- 
fect only the spectral region which is already contaminated 
by metal lines. Since P{k) is not affected by winds on large 
scales {k < 0.1 s km~^), the estimation of cosmological pa- 
rameters using the Lya forest flux power spectrum should 
be insensitive to galactic winds. 

(iii) The Lya flux transmissivity is enhanced near wind- 
blowing galaxies, independent of whether the wind is creat- 
ing a bubble of hot gas or a cooled cavity surrounded by a 
shell. This is consistent with the observations of AOS, who 
find that up to a third of all the galaxies in their sample with 
impact parameter smaller than lh~^ Mpc have transmitted 
fluxes up to 90%. The mean Lya transmissivity decreases 
with decreasing impact parameter, in agreement with pre- 
vious numerical models and the A05 data. 

(iv) To search for wind signatures in quasar spectra, we 
have devised a new method which measures the abundance 
of saturated regions as a function of their width. We find 
that absorption due to wind shells increases the expected 
number of saturated regions with _Ds < 1 A by a factor of 
two or more at 2; = 3 or 2, forcing the distribution to behave 
as a power law down to smaller widths than in the "no wind" 
model. In fact, the "no wind" model shows a turn-over at 
Ds ~ 1 A at 2 = 3 and at Db ~ 0.6 A at z = 2, while the 
wind models have such a point at Ds ~ 0.5 and Da ~ 0.2 
A respectively. Although not a huge effect, this is certainly 
large enough to be detectable. Since we consider saturated 
regions widths rather than line column densities, this mea- 
sure has the advantage of being insensitive to unsaturated 
metal lines blended with the Lya forest. 

(v) A similar analysis aimed to estimating the abundance 
of gaps in the absorption pattern of the Lya forest shows a 
weaker but still detectable effect from wind bubbles. We 
have calculated the cumulative distribution of gaps as a 
function of the gap width and we find that winds contribute 



as much as 50% more gaps with 10 < -Dg < 20 A in the eOl 
model. We suggest that this effect may be used to test the 
presence of hot winds in the Lya forest: if the cumulative 
distribution of gaps closely follows a power law, then it is 
likely that no winds are present, while if a bump appears 
for 10 < Dg < 20 A, then, according to our model, the Lya 
forest does contain wind signatures. 

(vi) The two new methods we present in Section|7|do not 
tell us directly where wind signatures are, but only whether 
they are present or not. However, this is already a significant 
improvement with respect to previous works. Finally, the de- 
tection of a large bump in the mentioned range may give us 
an indication that the powerful winds with low entrainment 
fractions of the eOl model are more likely to occur in the 
Universe than the weaker and more mass loaded winds in 
the e03 model. 
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APPENDIX A: 
STATE 



HYDROGEN lONISATION 



To calculate how much radiation of a desired wavelenght is 
absorbed by a gas cloud, we need to know the ionisation 
state of the gas, which determine its optical depth. To cal- 
culate the fraction of neutral hydrogen Xu I in a gas, we 
solve the equation for ionisation evolution of hydrogen. In 
the limit of i onisation equilibrium and for a highly ionised 
gas, we find JCeJl99l pTheuns et al.lll998^ : 



Xn I = 



r^H I + ToH iric 



(Al) 



where is the electron number density, an 11 the recom- 
bination rate of hydrogen, F^yH i the photoionisation rate 
of hydrogen and Fch i the collisional ionisation rate. The 
recombination rate of hydrogen is a function of the temper- 
ature of the gas and, de fining Tn = r/(10" K), is given by 
jHaardt fc Madaulll99el) : 



Oh II 



1.58- 10-"T4-'''^\ r<10*K 
1.58- 10-^^Tr°-5^-''-i'°'5'^^ T > 10'' K 



(A2) 

Hydrogen is ionised through two different mechanisms: 
collisional ionisation and photoionisation. Photoionisation is 
the dominant process in the low density gas with T < 10^ 
K. The photoionisation rate depends on the flux spectrum 
of the ionising UV background radiation J(i^, 2): 



F^H 1(2) 



4:TtJ{v, z)aniv) 
hv 



Av, 



(A3) 



where is the photoionisation cross-section and the 



Galactic Winds 13 



ionising threshold frequency of hydrogen. We use the ion- 
ising UV background J(i', z) and t he ph otoionisation rate 
r^H i{z) of'ifaard t. Madau fc R.ee3 l|l999h. At temperatures 
higher than about 10^ K, ionisation by collisional excita- 
tion becomes the dominant mechanism. In this regime, th e 
collisional ionisation rate is JCenlll99llTheuns et al]ll99d) : 



r,Hi = 1.17 X 10- 



-157809. l/T 



3 -1 

cm s 



(A4) 



APPENDIX B: OPTICAL DEPTH AND FLUX 
SPECTRUM 

The optical depth of the gas is calculated by integrating 
the relevant quantities (density, temperature and velocity) 
along a given LOS. We identify a sightline through the high 
resolution sphere and we divide its total length L into A'' 
pixels of equal width A = L/N. We choose A'' = 3000. The 
density and the density-weighted temperature and velocity 
for each bin j at position x{j) are com puted by integratin g 
the relevant quantities along the LOS iTheuns et alJll99dl : 



(pT)x(j)=a^^X(i)>Vijr(i 



(Bl) 



(B2) 



where nx{j) is the numer density, Tx{j) the temperature in 
pixel j and c the speed of light. The Lya cross section for H 
I is 



= 4.45 • 10"^** cm^ (B8) 



with Ao = 1215.6 A the rest wavelength of the transition 
and ctt = 6.625 • lO"'^^ cm^ the Thomson cross section. / — 
0.41615 is the oscillator strength and measures the quantum 
mechanical departure from the classical harmonic oscillator. 
We normalise our spectra by rescaling the mean flux 



^) to the mean flux measured in real spectra, e 



th I is the H i opacity along our simulated LOS. The effective 
optical depth t^s derived from observations is calculated as 
a function of redshift by using the relati on b etween redshift 
and optical depth given bv lKim et al.l ll2004l) in Subsection 
3.3 of their paper. The effective optical depth Tass in our 
spectra is: 



ross = -ln(e-"Hi). 



(B9) 



With an iterative procedure we rescale Toss to TcS until 
the difference between the two is no larger than a few per- 
cent, that is, typically, |ross — Tefi| < 0.005. This procedure 
allows us to make the artificial spectra directly comparable 
to the observed ones in each given redshift range. In addi- 
tion, it somewhat washes out the dependence of the optical 
depth on the normalisation coefficient of the UV ionising 
background radiation J21, which determines the ionisation 
fraction of hydrogen and therefore the absorption of the UV 
radiation. 



where the sum is over all the particles intersecting the LOS. 
a is the scale factor and X{i) is the abundance of the species 
X for particle i, assuming ionisation equilibrium. For Lya 
absorption spectra, X is the neutral hydrogen H i. The nor- 
malised SPH kernel is: 

mWinJhi) 



with W{rij/hi) given by 

„ ( 1-61 



W {v/h) = 



r/i3 



2(1- 
0, 



7^)3 + 6(^)^ 



(B4) 



i < ^ ^ 1> 
r > 1. 



(B5) 

as used in GADGET I. m is the mass of the dark matter 
particles, r^j is the distance between particle i and particle 
j and hi the smoothing length of the dark matter particle. 

In redshift space, a pixel at velocity v{k) suffers absorp- 
tion from a pixel at velocity v{j) by an amount e"^*-*^ , where 
r(fc) is the contribution to the optical depth of bin k given 
by bin j: 



r{k) 



nx{j)aA ■ exp 



v{k) - v{j) 



(B6) 

The Doppler width Vx{j) of the species X with mass 
mx determines the thermal broadening of the absorption 
lines. The Doppler width Vx(j) in pixel j is given by: 

2fcBrx(j) 



vm = 



mx 



(B7) 
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